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SUMMARY 


Expressions based on linearized supersonic-flow theory are derived 
for the perturbation velocity potential in space due to wing thickness 
for rectangular wings with biconvex airfoil sections and for arrow, 
delta, and quadrilateral wings with wedge-type airfoil sections. The 
complete range of supersonic speeds is considered subject to a minor 
aspect-ratio — Mach number restriction for the rectangular plan form and 
to the condition that the trailing edge is supersonic for the sweptback 
wings. The formulas presented can be utilized in determining the 
induced-flow characteristics at any point in the field and are readily 
adaptable for either numerical computation or analytical determination 
of any velocity components desired. 


INTRODUCTION 


The increasing use of auxiliary bodies such as stores and missiles 
on current aircraft has emphasized the need for reliable methods and 
procedures for predicting the load distribution, forces, and moments 
generated by such bodies and also the aerodynamic effects induced on 
these bodies by neighboring airplane components. The total loading 
and forces acting on stores, missiles, pylons, and so forth are required 
in order to design supporting structures, to predict the performance 
and stability characteristics of aircraft, to determine the jettison 
characteristics of stores, and to compute the trajectories of missiles 
when fired from parent aircraft. 

Although considerable effort has been expended in the past in cal- 
culating the load distribution, forces, and moments acting on isolated 
aircraft components for the supersonic speed regime (e.g., refs. 1 to 3 
and reports cited therein), somewhat less attention has been focused 
on the equally important problem of determining the aerodynamic effects 
induced by one body on neighboring airplane components. (See refs. 4 
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to 6 and reports mentioned therein.) The problem of determining the 
flow fields arising from airplane components undergoing various motions 
becomes a prime consideration in calculating induced aerodynamic effects 
Some recent work pertaining to the calculation of flow fields at super- 
sonic speeds may be found in references 7 bo 9- In general, the avail- 
able literature covers in adequate detail the flow fields arising from 
flat wings at an angle of attack and flat vings undergoing steady 
rolling or steady pitching motions. 

The present paper concerns itself with the linearized- theory evalu- 
ation of the effect of wing thickness and thickness distribution on the 
flow fields generated by thin wings at zero angle of attack. The pur- 
pose is to present closed-form expressions for the perturbation velocity 
potential which in turn may be utilized to calculate by either numerical 
or analytical procedures the desired flow-field velocities and angular- 
ities. Wings of rectangular plan form with symmetrical biconvex profile 
and of delta and modified-delta plan forms with wedge-type airfoil sec- 
tions are treated in detail. The complete range of supersonic speeds 
is considered subject to a minor aspect-ratio — Mach number restriction 
for the rectangular wing and to the condition that the trailing edge is 
supersonic for the sweptback wings. 


SYMBOLS 


£,y,z 




6**1 

V 


M 


P 

0 

0i>0ii> • • - 0ix 


rectangular coordinates cf field points 
nondimens ional rectangular coordinates defined as 
57? and 57? respectively 

z-coordinate defining airfoil section 


rectangular coordinates cf source points 
free-stream or flight velocity 
Mach number 

Mach number parameter, \/m 2 - 1 
perturbation velocity potential 

specific evaluations of pertubation velocity potential 


A 


sweepback angle of leading edge 
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m 

6 

c ] 

b 

S 

A 

A 

t 


sweepback-Mach number parameter, 0 cot A 

inclination of wing trailing edge, measured relative 
to root-chord extended (see fig. 2) 

root chord of wing 

wing span 

area of wing plan form 
wing aspect ratio, b^/s 

aspect-ratio — Mach number parameter, A0 
maximum thickness of local airfoil section 


c local wing chord 

A slope of wing surface, measured in stream direction 


F 1 * F 2 * * * - f 6 
Gf j * • • G i4 


Ha^ * * ,r lf 


functions that are given in appendix A for purposes 
of evaluating perturbation velocity potential for 
sweptback wing cases 

functions that are given in appendix B for purposes 
of evaluating perturbation velocity potential for 
rectangular wing cases 

specific values of variable t ) corresponding to 
various locations where Mach forecone from field 
point (x,y,z) intersects leading and trailing edges 
of wing (see fig. 3) 


ANALYSIS 


The analysis is based on supersonic thin-airfoil theory and on the 
assumptions of small disturbances and a constant velocity of sound 
throughout the fluid. These assumptions lead to the linearized equa- 
tion for the perturbation velocity potential 0: 


(1 - M2)0^ + 0yy + 0 ZZ = 0 


(1) 
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where M is the Mach number of the flow, and the derivatives are taken 
with respect to the variables x, y, and z of the rectangular coor- 
dinate system. The general expression for the linearized perturbation 
velocity potential in space due to a distribution of source and sink 
singularities in the z = 0 plane is (see refs. 10 and ll) 


0U,y,z) = - J f f — - - -ll = . (2) 

JJ \|(x - |) 2 - (m 2 - l)[(y - t !) 2 + z2] 

where x, y, and z are the rectangular coordinates of the field point 
at which the potential is desired, and £ and r\ are the rectangular 
coordinates (analogous to x and y) of the singularities. The func- 
tion A(|,tj) represents the particular dir.tr ibution of singularities 
for the wing under consideration and is thus, of course, dependent upon 
the boundary conditions imposed. For wing- thickness distributions that 
are amenable to thin- airfoil-theory calculations, the source- sink dis- 
tribution function is related to the particular thickness distribution 
involved and is given as 


Z B“ U 

The integration indicated in equation (2) is performed over the region R 
that is enclosed by the traces in the z = 0 plane of the Mach forecone 
emanating from the point (x,y,z) and by the- wing plan-form boundaries. 
(See fig. 1.) For purposes of convenience, the Mach number parameter 

3 = ^ - 1 rather than M itself will le used in the expressions to 

be developed; equation (2) may then be rewritten in the more familiar 
form 


S? z b (5,T| ‘ 


( 3 ) 


0(x,y,z) 


V 

Jt 




dT i 

2 - 0^(y - t]) 


2 


J3 2 z 2 


W 


In order to obtain closed-form express ions for the potential func- 
tion 0(x,y,z), it is necessary to define the slope function A(|,t]) 
for the particular wing under consideration (see eq. (3)) and then to 
integrate over the appropriate region R. The present paper considers 



5 


delta and nodif ied-delta wings with wedge airfoil sections and rectan- 
gular wings with symmetrical biconvex airfoil sections. (See fig. 2.) 
For the sweptback wings, the equation for the upper surface is given by 


z 



(5) 


where x is measured positive rearward from the y-axis. (See fig. 3*) 
The slope function A(£,t|) is, of course, constant for the sweptback 
wings and is given by 


A 


t 

2c 


( 6 ) 


For the rectangular wings, the equation for the upper surface may 
readily be derived as 


z = 2t -(l - - 
c \ c 

where x is measured positive rearward from the leading edge. The 
slope function A(|,t)), which for the rectangular wings is independent 
of r|, is found to be 




( 8 ) 


In order to determine the appropriate regions of integration, the 
sweptback wings have been subdivided into nine cases (see fig. 4), each 
case being characterized by distinct Mach trace — plan-form-boundary 
intersections. In an analogous manner, the rectangular wings have been 
subdivided into eight cases (see fig. 5). Thus, if the perturbation 
velocity potential is desired at an arbitrary point (x,y,z), the Mach 
trace must be determined and then the appropriate integration performed 
in accordance with equation (4). The sketches given in figures 4 and 5 
cover all possible cases that are required to determine completely the 
potential function for the wings under consideration. Note that the 
only restriction applicable to the sweptback wings is that the trailing 
edge is supersonic; for the rectangular wings, the side edges may not 
interact (expressed mathematically by the condition A3 ^ l). (See 
fig. 2.) It is of interest to point out that cases VII, VIII, and IX 
for the sweptback wings can occur only when the wing has a supersonic 
leading edge . 
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In order to define mathematically tie required regions of integra- 
tion, the T] coordinate of the point of intersection between the Mach 

forecone | = x - P \f(y - q) 2 + z 2 and each plan-form boundary must be 

obtained. The specific intersection points are indicated by the filled- 
in circles shown in figure and the associated q values are readily 
found by simple algebraic processes to be 


cot A 


o p l 

1 - p^cot^A L 


3 2 y cot A + p \J(x cot A + y) 2 + z 2 (l - 3 2 cot 2 A) 


(9) 




cot A 


p P 

1 - 3 cot A L 


3 2 y cot A - 3 \l ( x cct A - y) 2 + z 2 (l - 3 2 cot 2 A) 


( 10 ) 


1c 


tan & 


3 2 tan 2 6 - 1 


2 

x - c r + 3 2 y tan 5 - 3y [(x - c r )tan 5 + yj + z 2 (l - 3 2 tan 2 5) 


( 11 ) 


% = 


tan & 


3 2 tan 2 6 - 1 


tan & 


- x + tan 5 + 3^[(x - c r )tan & - yj + z 2 (l - 3 2 tan 2 S) 

(l 2 ) 


3 2 tan 2 & - lL 


’If 


cot A 


1 - p^cot^A 1 - 


c r -x + p 2 y tan 8 - Py[( x " c r )tan 8 - yj + z 2 (l - p^tan 2 8) 

(13) 

3 2 y cot A + 3 \|(x cct A - y) 2 + z 2 (l - 3 2 cot 2 A) 


( 1*0 


For the rectangular wing (A = 0° and 8 = 90°) > q e and are 

the same as r\ c and T] a , respectively, and considerable simplification 
results in equations ( 9 ) to (12). The simplified expressions for r) a 
to are as follows: 
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The integral expressions for the perturbation velocity potential 
corresponding to the nine cases for the sweptback wing (see fig. 4) and 
the eight cases for the rectangular wing (see fig. 5) may now be explic- 
itly determined. For purposes of convenience appendixes A and B present 
for the sweptback wings and rectangular wings, respectively, the various 
conditions associated with each potential function as well as the spe- 
cific formulas (F and G functions) required to define the potential. 


The following equations result for the nine sweptback-wing cases, 
where the slope A is constant and defined by equation (6): 


Case I: 


Case II: 


n 0l( x ,y,z) 

VA b/2 


* 0n(x,y,z) , w 

i? 1 + i*p 

VA b/2 1 d 


( 19 ) 

( 20 ) 


Case III: 


it 0ni( x ' y - >z ) 

VA b/2 


Ff + F2 + Fj 


(21) 


Case IV: 


jt_ 0iv( x ^y j z ) _ F + F 

VA b/2 1 


F, 


5 


( 22 ) 
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Case V: 


Case VI: 


Case VII: 


Case VIII: 


Case IX: 


_n_ 0v(*,y,z) 
VA b/2 


n_ 0vi( x ^y^ z ) 
VA b/2 


+ F2 + Fj + 


n_ 0vn(x,y,z) = 

VA b/2 5 


- f 5 + f 6 


jt 0 vm( x >y> z ) 

TO i/I F 6 ' F 5 


_n_ 0ix( x ^y^ z ) 
VA b/2 


(23) 

(24) 

(25) 

(26) 

(27) 


For the eight rectangular-wing cases, the following formulas 
result : 


Case I: 


n 0i( x iy>z) 

Kt) b/2 


- 5G]^ + 


(28) 


Case II: 


For y < — , 
J " 2 J 


For y > — , 


2VI 


It 01l( x ;y, z ) 

7t\ ~ V2 = Gi - G 3 + G 4 + G 5 


n 0n( x >y> z ) 


> (29) 


2 V 1 


t \ b/2 


— G2 - Gj + G4 


Case III: 


it 0 m( x >y> z ) 


Zv{±\ b /2 


= G]_ + G;_q - Gj + G^i_ + G^ 


(30) 
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Case IV : 


n 0rv( x >y; z ) 

2vf—) b / 2 


(3D 


Case V: 

5— - = G i + G io " g 3 + G n + g 5 " G 6 ' g 12 + g 8 “ g 13 

«(!) 


Case VI: 


Case VII: 
For 

For 

Case VIII: 


2 VI 


_n 0yi(x,y,z) 

~I\ D/2 ” 


= Gi + G]_q - g 3 + G 11 “ 2G4 


y i 


< £ 


2' 


Jt _ r n, n_2 

b y 2 G x - Gg - G3 + GQ + g 9 


y > 2 > 


2 v(i) 

it 0vn( x ^y^ z ) 


-d) 




= G2 ” Gy - Gj + Gg - G^ 


0vm( x ^y^ z ) 


K!' 


b/2 


= G]_ - G j - G4 


G lU 

(32) 

(33) 


(34) 


(35) 


DISCUSSION 


The formulas presented in the previous section enable the direct 
evaluation of the perturbation velocity potential in the vicinity of 
rectangular wings with symmetrical biconvex airfoil sections and the 
corresponding perturbation velocity potential for delta and modified- 
delta wings with wedge airfoil sections. For the sweptback-wing cases 
with airfoil sections of arbitrary thickness distribution, the well- 
known superposition procedure can be employed by using the given formulas 
with appropriate values for the parameters substituted therein. A graph- 
ical presentation of the mechanics involved in superposing the basic 
solution is given in figure 6. Note that each component part may be 
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directly evaluated from the expressions presented in equations (6) and 
(19) to ( 27 ) and those given in appendix A. 


In connection with the formulas applicable for the sweptback wings, 
it should be pointed out that for the special conditions of (3 cot A = 1 
(sonic leading edge) and p tan 5 = 00 (unsvept trailing edge) the func- 
tions of appendix A are simplified to consilerable extent by conventional 
procedures. The following relationships should be used in the event 
certain terms (for convenience denoted by I'J) take on imaginary values: 

For | N | < 1, 


and, for N > 1, 


cosh = i cos~^TJ 


cos'^N = -i cosh ^ 


cosh ^ = -i cos~% 
cos"^N - i cosh~^N J 


( 36 ) 


( 37 ) 


Care should be exercised in extracting only positive roots from radical 
terms. To minimize the possibility of errcrs in the final formulas, all 
derivations have been checked analytically by independent means and also 
checked numerically by graphical integration procedures. 


The formulas presented for the perturbation velocity potential may 
be differentiated by either numerical or analytical procedures to obtain 
the various velocity components which in turn will enable the direct 
evaluation of flow-field effects on neighboring airplane or missile 
components. An approach that has proved to be very efficient consists 
of computing the potential function and its variation along a given 
direction and then measuring the slopes either graphically or numerically 
to obtain the desired velocity component. 


CONCLUDING REMARKS 


Expressions based on linearized supersonic-flow theory are derived 
for the perturbation velocity potential in space due to wing thickness 
for rectangular wings with biconvex airfoil sections and for arrowhead, 
delta, and quadrilateral wings with wedge-type airfoil sections. The 
complete range of supersonic speed is considered subject to a minor 
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aspect -ratio— Mach number restriction for the rectangular plan form and 
to the condition that the trailing edge is supersonic for the sweptback 
wings. The formulas presented can be utilized in determining the induced- 
flow characteristics at any point in the field and are readily adaptable 
for either numerical computation or analytical determination of any 
velocity components desired. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., December 2b, 1958. 
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APPENDIX A 

MATHEMATICAL CONDITIONS AND FUNCTIONS PERTINENT TO EVALUATION 
OF VELOCITY POTENTIAL FOR SWEPTBACK-WING CASES 

In order to determine the specific sweptback-wing case applicable 
for a given point (x,y,z), the following mathematical conditions may be 
utilized: 

Case I: 

^a < 0 < % 

T) c and are imaginary 


Case II: 


Case III: 


Case IV : 


Case V: 



0 < < \ < 

0 < % < ^d < ;; 
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Case VI: 



Case VII: 

0 < n f < Tie < | < v 

Case VIII: 

0 < n f < n e < % < I 

Case IX: 

0 < nr < % < I 

Tj e and are Imaginary 

The various r\ functions referred to in the preceding conditions are 
defined by equations (9) to (l 4 ) in the text. 

The formulas for the perturbation velocity potential 0 are pre- 
sented for the various cases in equations (19) to (27) in the text. 

The functions F^_ to Fg are used therein for purposes of simplifica- 
tion and are defined as follows in terms of the nondimens ional space 
coordinates x, y, and z and the plan-form — Mach number parameters A 
and m (where A = Ap and m = p> cot A) : 




Ik 


t 2 = - 


M* - i) + yfj - M cos -i 


4(x - l)(A - 4m) + Tftny 




A 2 - 16 ^ + 8 Am - A 2 


l|fM* - 1) + y(I - 4m)] 2 - Z^m^A 2 - 16 ) + 8Am - A 2 


Q -1 rn£^/l6(x - l) 2 - A^y 2 + + - tan -i i^^6(x - l) g - ^(y 2 + z 2 ) 

(A - 4n )(y 2 + z 2 ) - 4my(x - l) 


(A - 4m)(^y 2 + z 2 j + Mmy(x - l) 
4m(x - 1) - y(A - 4m) x 


4 ( x - l) (/ - 4m) - A 2 my 


jjm 2 (^ - l6j + 8Am - A 2 A |[4m(x - l) - y(A - 4m)] ‘ - z 2 m 2 (A 2 - l6) + 8Am - A 2 ] 


(A2) 


F< - 


4mx - Ay , -1 4mx - A + Am 2 ( 1 - y) 
— cosh 


z\l(h m* - A) 2 - sVfa - y) 2 + i 2 ] 


S '/ 1 - m2 m||(innJJ - lyf + 7?Z 2 (l - m2) 5 [s 2 - y(l - ?)] + ^(l - ?) 

4m(x - l) - y(A - 4m) (A - 4m ) ( 4mx - A) + A 2 m^(l - y) _ 


ijm^A 2 - lfi) + 8Am - A 2 Am y [4m( x - 1) - y(A - 4m )_ - z 2 rr^A 2 - l6) + 8Am - A 2 ! 


z^(bi 

nx - A) 2 - A^m 2 

[a - y ) 2 + 2 ? _ 


(A - 4m) 

_z 2 - y(i - y)_ 

+ 4m(x - 1)(1 

- ?) 


(A3) 


Fji - 


4m(x - l) + y(A - 4m) CQG -1 
a2 (a 2 - l6^) + 8 Am - A 2 


(A - 4m)(4mx A) + A 2 m 2 (i + y) 


Am 


^ [4m(x - 1) + y(A - 4m) - z 2 jm 2 (A 2 - l6^ + 8Am - A 2 ^ 


S \ 

( 4mx - A) 2 - A 2 [m 2 (l + y) + z 2 ] 

4m; + Ay 1 4mx - A + Ain 2 (l + y) 

3 

i 

z 2 + y(l + y) 

+ 4m(x - l)(l 

+ y ) Al|. - m2 m|(4mx + Ay ) ^ + A 2 z 2 (l - m2) 


z tan 


^ z|(4mx - A) 2 - A^m 2 [( 1 + y) 2 + z 2 ] 
A z 2 + y(l + y) + 4mx(l + y ) 


(A4) 


4m(x - l) - y(A - 4m) 


^(a 2 - 16) 


+ 8 Am - A 


F 6 = feL=£- Z|« 

l A||m2 - 1 


(A5) 


(A6) 





15 


APPENDIX B 

MATHEMATICAL CONDITIONS AND FUNCTIONS PERTINENT 
TO EVALUATION OF VELOCITY POTENTIAL FOR 
RECTANGULAR-WING CASES 


In order to determine the specific re ct angular- ving case applicable 
for a given point (x,y, z) , the folloving mathematical conditions may be 
utilized: 


Case I: 


x < c + (3z 



% < % < | 


Case II: 


Case III: 


For 

For 


x < c + 0z, 
x > c + J3z, 


For x < c + 


For x > c + 


' I < la < ! < % 
‘ I < V < ! < 1c 


0z, 




b 

2 


% > 


b 

2 


Pz, 


< - 


b 

2 



Case IV: 


x > c + pz 
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Case V: 



Case VI: 

v 

^a < “ 2 < ^ < ^ < 2 < ^ 

Case VII: 

- | < % < < % 

Case VIII: 

- 1 < n a < n c < % < | < % 


The various r| functions referred to in the preceding conditions are 
defined by equations (15) to (l8) in the text. 

The formulas for the perturbation velocity potential 0 are pre- 
sented for the various cases in equations (28) to (35) in tl 16 text. The 
functions Gj to are used therein for purposes of simplifi cation 

and are defined as follows in terms of tie nondimensional space coordi- 
nates x, y, and z and the aspect-rat Lo — Mach number parameter A 
(where A = Ap): 
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G 2 = (l - 2x) I (l - y)cosh -1 — — — + 

Ap - l) 2 + z 2 

\ 

£ cos _ 1 z 2 (ux 2 - A^z 2 ) - (y - l) 2 (lx 2 + A^z 2 ) 
^4x 2 - T^z 2 ) (y - l) 2 + z 2 J 

x -1 2l 2 (y - l) 2 - ^x 2 + A 2 z 2 

— COS 

A 4x 2 - A 2 z 2 




G 5 = (1 - 2 x)%^ 
A 


(B5) 





(b6) 
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Gj = (l - 2x) < (l - y)cosh -1 — !===== = — + 

[ A^(y - l) 2 + Z2 

z i z 2 t-( s - l) 2 - A 2 z 2 j - (y - l) 2 [4(x - l) 2 + A 2 ^ 
— cos -J - 1= - — 

[4(x - l) 2 - A 2 Z 2 ] [(y - l) 2 + z 2 ] 

x - 1 _ L 2ff 2 (y - l) 2 - 4(x - l) 2 + I 2 z 2 
A 4(x - l) 2 - X 2 z 2 





2x)— 

A 


(B9) 


G 10 = (1 - 2x) J(l + y)cosh -1 — = - 

} Sp + l) 2 *• 2 2 

g z 2 (4x 2 - A 2 ^ 2 ) - (y + l) 2 (4c 2 + A 2 ^ 2 ) 

2 (4x 2 - ^z 2 )[(y + l) 2 + z 2 ] 

5 _i ^(y + l) 2 - 43c 2 + A 2 ? 2 
A 4x 2 - A 2 z 2 


(BIO) 
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J li 


, L± . ^[(y + l) 2 + 3 + U * Z - t£ sin-1 A <? i X L- (Bll) 

2 v L ® ^ T ^ 7 


-t- 2_2 
A z 


g 12 = (l _ 2x) <(l + y)cosh 


-1 2 (x - 1 ) 


Al(y + I)* + 


2 ^ *2 


z -1 
— cos x 

2 


1 z 2 [4(S - l) 2 - A 2 ^ 2 ] - (y + 1) 2 [,U - l) 2 + A^ 2 ] 
_4(x - l) 2 - S 2 ^ (y + l) 2 + z 2 


x - 1 2A 2 (y + l) 2 - 4(x - l) 2 + A 2 ^ 


COS 


fl 


4(3? - l) 2 - 7?z 2 


(B12) 


g 13 = 


L±l\Jk(x - l ) 2 - ^[(y + l ) 2 + z 2 ]” + 

4(3? - l) 2 - A 2 ^ 2 . A(y + 1) 

— — sm ± - 

2A j^4(x - l) 2 - ^Z 2 


(B13) 


G l 4 = (1 - 23c) (x - 1)— 
x A 


(Bl4) 
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and the plan form solid; note that the region R is defined by the boundary abed. 
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